home *** CD-ROM | disk | FTP | other *** search
/ Over 1,000 Windows 95 Programs / Over 1000 Windows 95 Programs (Microforum) (Disc 1).iso / 1249 / newton.t < prev    next >
Encoding:
Text File  |  1997-04-18  |  1.0 KB  |  55 lines
  1. %
  2. % "newton.t" uses Newton iteration to solve a 
  3. % simultaneous nonlinear set of equations:
  4. %
  5. %       f1(x,y) = x^3 - x^2 * y + y^2 - 3.526
  6. %       f2(x,y) = x^3 + y^3 - 14.911
  7. %
  8. %   Sample program for the T Interpreter by:
  9. %
  10. %   Stephen R. Schmitt
  11. %   962 Depot Road
  12. %   Boxborough, MA 01719
  13. %
  14.  
  15. const DIM : int := 2
  16.  
  17. program
  18.  
  19.     var p, q : rmatrix
  20.     var x : real := 1.0
  21.     var y : real := 2.0
  22.     var det, f1, f2 : real
  23.  
  24.     var i : int
  25.     
  26.     for i := 1 ... 5 do
  27.  
  28.         p[0,0] := 3 * x^2 - 2 * x * y
  29.         p[0,1] := 2 * y - x^2
  30.         p[1,0] := 3 * x^2
  31.         p[1,1] := 3 * y^2
  32.  
  33.         
  34.         det := invert( p, q, true ) 
  35.  
  36.         if det ~= 0.0 then
  37.  
  38.             f1 := x^3 - x^2 * y + y^2 - 3.526
  39.             f2 := x^3 + y^3 - 14.911
  40.  
  41.             x := x - q[0,0] * f1 - q[0,1] * f2
  42.             y := y - q[1,0] * f1 - q[1,1] * f2
  43.  
  44.             put "x = ", x:12:8 , ", ", "y = ", y:12:8
  45.  
  46.         else
  47.  
  48.             put "singular"
  49.             exit
  50.  
  51.         end if
  52.  
  53.     end for
  54.  
  55. end program